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We present specific-heat measurements in the vicinity of the antiferromagnetic phase transition on single crystals of 
the alloy Yb(Rhi_2;Co K )2Si2 for x < 0.38. This study was motivated by the violation of critical universality in the 
undoped YbRli2Si2 (Krellner et al., Phys. Rev. Lett. 102, 196402) where we have found a large critical exponent 
a = 0.38. For Co-doped samples we observe a drastic change of the critical fluctuations resulting in a negative a, 
explainable within the universality classes of phase transitions. The development of a under chemical pressure gives 
strong indication that the violation of critical universality in YbRli2Si2 is due to the nearby quantum critical point. 
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1 Introduction In recent years, there have been con- 
siderable efforts to explore the physics of quantum phase 
transitions. These transitions are driven by quantum fluc- 
tuations in contrast to classical phase transitions which are 
actuated by temperature. In this context the heavy fermion 
system YbRli2Si2 was intensively studied, because it is a 
clean and stoichiometric metal situated on the magnetic 
side of, but very close to, a quantum critical point (QCP) 
which can be crossed by applying a tiny magnetic field. 
Therefore, this system presents both an antiferromagnetic 
(AFM) phase transition driven by thermal fluctuations as 
well as pronounced quantum fluctuations (for an experi- 
mental and theoretical review see e.g., Ref. HI and J2)). 
The recent discoveries of an additional energy scale van- 
ishing at the QCP which does neither correspond to the 
Neel temperature nor to the upper boundary of the Fermi- 
liquid region J5) and a large critical exponent a = 0.38 
at the AFM phase transition observed in low-temperature 
specific-heat measurements on a single crystal of supe- 
rior quality (4) have once again boosted the interest in 
YbRli2Si2. The latter observation triggered strong theoreti- 
cal effort to explain the violation of critical universality due 
to a (quantum) tricritical point 05] • In this scenario, Mi- 
sawa et al. J2) proposed the existence of a tricritical point 
for YbRli2Si2 under pressure and magnetic field at finite 



temperatures. Experimentally, this part of the phase dia- 
gram is easiest to explore using chemical pressure as will 
be discussed below. 

The magnetic ordering of YbRli2Si2 (T/v = 72 mK) is sta- 
bilized by applying pressure as expected for Yb-Kondo lat- 
tice compounds |6|. The complementary method of substi- 
tuting smaller isoelectronic Co for Rh results in chemical 
pressure allowing a detailed investigation of the magnetic 
phase diagram and the physical behavior of the stabilized 
AFM ordered state. Therefore, a thorough understanding 
of the physical properties of Yb(Rhi_ :r Co 2 ;)2Si2 is of great 
interest in order to get further insight into the phenomena 
at the QCP in YbRli2Si2. Very recently, it was shown that 
for x = 0.07 the signature of the Kondo breakdown is lo- 
cated within the magnetically ordered phase leading to a 
detaching of the AFM QCP from the Fermi-surface recon- 
struction [7j. In this contribution, we address the question 
whether the critical fluctuations around Tjv remain anoma- 
lous when moving away from the Kondo-breakdown QCP, 
using positive chemical pressure. To this end, we present 
detailed specific-heat measurements around TV for x = 
0.12, 0.27, and 0.38. 



2 Experimental Single crystals of the alloy series 
Yb(Rhi_ a; Co a ;)2Si2 were grown from In flux, analogous 
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Figure 1 4/ increment to the specific heat plotted as 
C 4 ' /T(T) on a logarithmic T scale for three single crys- 
tals of the series Yb(Rhi_. r Co :r )2Si2 with x = (circles, 
0), 0.12 (triangles), and 0.38 (squares), respectively. 



to the stoichiometric samples of superior quality |4]. The 
random substitution of Rh with Co leads to larger disor- 
der compared to YbRli2Si2; however, a comparison with 
pressure experiments on undoped YbRli2Si2 prove that 
Co-doping acts mainly as chemical pressure J7). A thor- 
ough investigation of the complete doping series, including 
x-ray diffraction, magnetic susceptibility, electrical resis- 
tivity, and specific-heat measurements was performed and 
will be published separately (8). The Co-content was ac- 
curately determined using energy dispersive x-ray spectra 
of the polished single crystals. This real Co-content will 
be used for x throughout the manuscript. The specific heat 
for the doped samples was determined with a commercial 
(Quantum Design) physical property measurement system 
(PPMS) equipped with an 3 He-insert, using a standard 
heat-pulse relaxation technique. The 4/ contribution to 
the specific heat, C 4 -', was obtained by subtracting the 
non-magnetic one, Cl u , from the measured specific heat, 
C m eas- C*Lu was determined by measuring the specific 
heat of the non-magnetic reference sample LuRli2Si2 be- 
low 10 K 1H. Since C Lu at 1 K contributes only to 1 % of 
Cmeas, the synthesis and measurements of the appropriate 
reference systems Lu(Rhi_ 2; Co a ;)2Si2 was not necessary. 

3 Results In Fig.Q] the 4/ increment to the specific 
heat is plotted as C if /T on a logarithmic temperature 
scale for three selected Co-concentrations, x = 0, 0.12, 
and 0.38. The anomalies due to the onset of the magnetic 
order are clearly visible and the type of the anomalies gives 
already a clue to the development of the critical fluctua- 
tions with increasing x. Going from x = to x = 0.12, we 
observe a drastic change of the peak form at Tjy. For the 
former we observe a very sharp A-type peak which leads to 
a large critical exponent |4|, whereas the latter presents the 
case of a rather broad mean-field-like anomaly, where it is 
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Figure 2 Specific heat vs reduced temperature close to 
Tjv of Yb(Rho.73Coo.27)2Si2- The data below and above 
T N can be best fitted with a 27 = -0.06 ± 0.10 (solid 



lines). Dotted lines represent the fit including a Gaussian 
distribution of T N with 6T N /T N = 6.5 • 10~ 3 with other- 
wise identical fit parameters. 



not possible to extract the critical exponent. For x = 0.38 
we again observe a sharp A-type anomaly which we can an- 
alyze in terms of critical fluctuations. This development to- 
wards a sharper anomaly with increasing x gives evidence 
that the broad anomaly for x — 0.12 is not caused by dis- 
order due to the random distribution of the Co-atoms on 
the Rh side, as the residual resistivity increases together 
with the doping level from x = 0.12 to x = 0.38. It is 
important to note, that for 0.07 < x < 0.18 we observe 
two subsequent magnetic phase transitions ifTOl . whereas 
for 0.27 < x < 0.38 we see only one magnetic transi- 
tion {8) indicating that the magnetic ordering vector might 
change with increasing x. 

Next, we focus on the analysis of the critical exponent 
for the two Co-concentrations, x = 0.27 and 0.38. For 
x = 0.27 the peak at T/v looks similar to what is presented 
for the x = 0.38 sample (8). To extract the critical expo- 
nent we follow the procedure described in detail in Ref . J4] . 
The usual fit function is applied to describe the critical be- 
havior 

A ± 



C 



: (t) = —\t\ 
a 



b + Et. 



(1) 



with the reduced temperature t = (T — Tjv)/Tjv; +(— ) 
refers to t > (t < 0), respectively ifTTIl . The background 
contribution is approximated by a linear t dependence 
(b + Et) close to T N 0). The best fit for x = 0.27 reveals 
a negative critical exponent a 27 = —0.06 ± 0.10 (Tab.[T]i 
and is shown together with the experimental data in Fig. [2] 
To satisfactorily fit the data points for \t\ < 0.01 we have to 
use a Gaussian distribution of T/v with STn = 8.4 mK to 
account for the rounding effects, most probably due to the 
effect of the high doping level. This smearing of Tjv leads 
to a relatively small temperature range, 0.01 < \t\ < 0.1, 
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Figure 3 Specific heat vs reduced temperature close to 
Tpf of Yb(Rho.62Coo.38)2Si2- The data below and above 
T N can be best fitted with a = -0.12 ± 0.10 (solid lines). 
Dotted lines represent the fit including a Gaussian distri- 
bution of Tjv with 5Tn /Tn = 5 • 10~ 3 with otherwise 
identical fit parameters. 

Table 1 Parameters obtained from the fits of the specific- 
heat around Tjv for x = (4)> x = 0.27, and x = 0.38. 



X 


Tjv 


STn /Tn 


A+/A- 


a 




(K) 


(io- 3 ) 









0.072 


0.3 


0.6(1) 


+0.38(3) 


0.27 


1.298 


6.5 


1.4(5) 


-0.06(10) 


0.38 


1.223 


5 


2.2(5) 


-0.12(10) 



which determines the critical exponent, severely compli- 
cating the analysis and impeding a more accurate deter- 
mination of a. The same analysis was carried out for the 
measurement of the x = 0.38 sample and is shown in a 
similar way in Fig. [5] The curves exhibit a comparable 
overall i-dependence like the x = 0.27 data, giving again 
a negative critical exponent a 38 = —0.12 ± 0.10. Here, 
we have to use a slightly smaller distribution of Tn with 
8Tn/Tn = 5 • 10~ 3 consistent with a lower residual resis- 
tivity in the x = 0.38 crystal (H. 

In Tab. [TJ we summarize the parameters obtained from 
the analysis of the critical exponent for the doped sam- 
ples presented here and compare them with the ones ob- 
tained for the undoped compound [4]. It is obvious that 
there is a drastic change of the critical exponent going 
from a = 0.38(3) for x = to a = -0.12(10) for 
x = 0.38. The latter value can be explained in terms of 
the classical universality classes in the theory of phase 
transitions, for which -0.133(5) < a < +0.110(1) gen- 
erally holds true. However, the low experimental accuracy 
of the determined exponent is not sufficient to distinguish 
between the two applicable symmetry classes ifTTl . namely 
the 3D-Heisenberg model [a^o.H = —0.133(5)] or the 
3D-XY model [a^ D:X Y = -0.015(1)]. More important 



is the clear development of a from an unconventional 
value to a conventional one when moving away from the 
Kondo-breakdown QCP by chemical pressure which gives 
a first indication that the violation of critical universality 
in YbRli2Si2 is due to this nearby unconventional QCP 
which may substantially influence the spatial fluctuations 
of the classical order parameter. 

Furthermore, our result of a conventional critical exponent 
for YbRli2Si2 under chemical pressure is in contradic- 
tion to what is expected in the quantum-tricritical-point 
scenario Q, where Misawa et al. suggest a stabilized 
tricritical point at finite temperatures in the pressurized 
YbRh2Si2- In this case one would expect that the critical 
exponent becomes even larger for YbRli2Si2 under chemi- 
cal pressure, as the theoretical value at a classical tricritical 
point is a = 0.5, just opposite to what we have observed. 

4 Conclusions In conclusion, we have presented 
specific-heat measurements of Co-doped YbRli2Si2 in 
the vicinity of Tn and studied the development of the 
critical exponent. We found that the critical fluctuations 
change drastically with doping. The broad mean-field- 
type anomaly at Tn for 0.07 < x < 0.18 prevents the 
determination of the critical exponent, but for x = 0.27 
and 0.38 we found a conventional critical exponent of 
a = —0.06(10) and —0.12(10), respectively. The develop- 
ment of a from an unconventional value to a conventional 
one when moving away from the local QCP by chemical 
pressure gives strong indication that the violation of critical 
universality in YbRli2Si2 is due to this nearby QCP. 
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